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A thickening at the outer edge of the accretion disk is usuahy invoked to explain the 
dips in the light curves of cataclysmic variables with stationary disks at phases ~ 0.7. The 
non-collisional interaction between the stream and the disk in the stationary solution raises 
' the question of why matter appears at a considerable height above the accretion disk in 

such systems. Our three-dimensional numerical modeling demonstrates that a thickening 



10 . of the halo above the disk can appear even in the absence of a direct collision between 

^2 I the stream and the disk. In the gas-dynamical flow pattern described with the "hot- 

f-*) ■ line" model, a considerable fraction of the matter is accelerated in the vertical direction 

10 . during the flow's interaction with the circumdisk halo. The vertical motion of the gas due 

I to the presence of the z component of the velocity leads to a gradual thickening of the 

i-C ■ circumdisk halo. The computations reveal the strongest thickening of the halo above the 
outer edge of the disk at phases ~ 0.7, in agreement with observations for stationary-disk 

O I cataclysmic variables. This supports the hot-line model suggested earlier as a description 

' of the pattern of the matter flows in semi-detached binaries, and presents new possibilities 

■ for interpreting the light curves of such systems. 

> ' 

^ 1 Introduction 

_ _ I 

Observations of low-mass X-ray binaries (LMXBs) have revealed dips in the X-ray light 
curves for several systems. Explanations of these dips have hypothesized the presence of a 
thickening at the accretion disk's outer edge at phase ~ 0.8, which corresponds to the position 
of this feature in the light curves ([1, 2] and references therein). The presence of matter 
surrounding the X-ray source at a considerable height above the system's orbital plane, along 
with the matter's uneven distribution in azimuth, can be explained in terms of either the 
companion's gravitational action on the accretion disk or the interaction between the stream 
of matter from the inner Lagrangian point (Li) and the disk. The coincidence of the phase 
of the observed dips with the assumed position where the stream from Li approaches the 
outer edge of the disk has tended to focus study on this particular region. Beginning with 
[3-5], numerical studies of variations of the disk's scale height due to its interaction with the 
stream were initiated. Currently, the idea that the stream ricochets off the disk's outer edge 
is considered the best way to explain the presence of matter at heights considerably exceeding 
the disk thickness. The possibility that the stream flows around the edge of the disk was first 
discussed in [6-9]; in 1996, the gas-dynamical computations of Armitage and Livio [10] showed 
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Figure 1: Light curves of the UX UMa (left) and RW Tri (right) systems in the ultraviolet (1415 ^ 
I433A), according to [16]. 



that a considerable fraction of the stream's matter can ricochet off the edge of the accretion 
disk. In their model, after colliding with the disk edge, some of the matter in the stream rises 
to a considerable height (compared to the disk thickness), forming a stream towards the inner 
parts of the disk. The computations of [10] demonstrated that this stream of matter could 
explain the presence of the dips observed in LMXB light curves. 

It is easiest to identify dips due to the presence of matter located high above the disk 
in LMXBs, since they contain a very compact source at the center of the disk. However, 
similar light-curve features in various wavelength ranges have also been recorded for a number 
of cataclysmic binaries in outburst, such as U Gem [11, 12], OY Car [13, 14], and Z Cha 
[15]. Further observations showed that lightcurve dips can also appear when a system is in a 
stationary state. Studies of the ultraviolet hght curves of the eclipsing nova-like cataclysmic 
binaries UX UMa and RW Tri [16] confirmed this result, and suggested that this phenomenon 
was universal in semi-detached binaries with accretion disks. As an example. Fig. ^ presents 
the observed light curves for UX UMa (left panel) and RW Tri (right panel) from [16]. It is 
interesting that, in contrast to the cataclysmic systems, systems with stationary disks display 
pre-eclipse dips at much earlier phases^, 0.6 -7- 0.7 [16, 17] (Fig. [T)). 

This raises the question of what leads to the presence of matter at considerable heights 
above the accretion disk in the case of a stationary interaction between the stream and the 
disk. Gasdynamical studies of the terminal flow pattern in semi-detached binaries demonstrate 
that the interaction between the stream and the disk is coUisionless in this case [18-24]. In 
contrast to the hot-spot model, which assumes that the stream impacts the edge of the accretion 
disk, the stream interacts with gas of the circumdisk halo in the stationary case, forming an 
extended region of enhanced energy release, or so called "hot line". This means that there will 
be no ricochet of the stream in the steady-state mode, and hence this model cannot explain 
the dips in the light curves of binaries with stationary disks. 

Our aim here is to study possible ways of thickening the halo above a disk, which give rise to 
eclipses of the central source and the presence of dips in the light curves during the stationary 
flow of matter in semidetached binaries. Our three-dimensional modeling of the flow structure 
demonstrates that, during the stream's interaction with the circumdisk halo, a considerable 
fraction of the matter is accelerated in the vertical direction and rises to heights appreciably 
exceeding the disk scale height. The strongest thickening of the halo above the outer edge of 
the disk occurs at phase ~ 0.7, in agreement with observations of systems with a stationary 
disks. This provides new possibilities for interpreting the light curves of such systems. 

"'^As usual in the analysis of observational data, the phase angle, (j>, is measured from the line connecting the 
centers of the two stars, with the phase increasing in the direction opposite to the system's rotation. 
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2 The Model 



We used the model of [22] in our numerical study of the gas dynamics of the matter flows 
in a semidetached binary. We described the flow pattern using a three-dimensional system of 
gravitational gasdynamical equations including the effects of radiative heating and cooling of 
the gas: 

^ + divpv = 0, 

^ff- + div(pv (8) v) + grad P ^ -p grad $ - 2[f2 x v]p , 
< (1) 

dp{e +^MV2) ^ + + |^|2/2) = 

-pv grad $ + p^m-^ {V{T, T^a) - A{T)) . 

Here p is the density, v = (u, v, w) the velocity vector, P the pressure, e the internal energy, 
$ the Roche potential, rrip the proton mass, and r(T,Tyjd) and A(T) the radiative heating 
and cooling functions, respectively. The system of gas-dynamical equations was closed with the 
equation of state for an ideal gas, P = {■j — l)ps, where 7 is the adiabatic index. The parameter 
was taken to be 5/3. 

We solved this system of equations using the Roe-Osher technique [21, 25, 26] adapted for 
multiprocessor computers. The modeling was carried out in a noninertial frame of reference 
rotating with the binary, in Cartesian coordinates, on a rectangular three-dimensional grid. 
Since the problem is symmetric about the equatorial plane, we modeled only half of the space 
occupied by the disk. The size of the modeled region, 1.12^4 x 1.14^4 x 0.17^4 (where A is the 
distance between the system's components), was chosen to completely include the disk and the 
stream of matter from the point Li. The computational grid had 121 x 121 x 62 cells distributed 
among 81 processors forming a 9 x 9 two-dimensional array. The grid was made denser in the 
region of interaction between the stream and the disk in order to improve the accuracy of the 
solution in this region. The grid was also denser near the equatorial plane, providing good 
resolution of the disk's vertical structure. 

The solution obtained for the model without cooling was used for the initial conditions [27]. 
Before the solution converged, the model with cooling was computed for ^ 5 orbital periods 
of the system as a whole. The total computation time was ~ 1000 hours at the MVS lOOOA 
computer of the Joint Supercomputer Center (JSCC). 

A free boundary condition with constant density p^ = 10"^^^,^, a temperature of 13600° K, 
and zero velocity was imposed at the outer boundaries of the computational domain (except 
near Li ), where Li is the matter density at Li. The accretor was taken to be a sphere 
of radius 10~^A. All the matter entering the accretor cells was taken to fall onto the star. 
The stream was specifled as a boundary condition: matter with temperature 5800° K, density 
Pl^ — 1.6 X 10~^g/cm^ , and velocity along the x axis — d.Skm/s was injected into a region 
with radius 0.0 14 A around Li. 

We considered a semi-detached binary consisting of a donor with mass M2 that fllls its 
Roche lobe and an accretor with mass Mi. The following parameters were adopted for the 
system: Mi = I.O2M0, M2 = O.5M0 , and A = 1.42i?0 , corresponding to Port = 3.79'*. The 
disk in the model had a temperature of 13600°ir. For the specifled rate of matter entering the 
system, the corresponding accretion rate in the model was ~ lO~^°M0/yr. 
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Figure 2: Three-dimensional surface of constant density p = 5 x 10 ^^gjcrrc' . 

3 Results of the Computations 

We described the morphology of the matter flows in a semi-detached binary with a sta- 
tionary, cool (T = 13600°ii') disk in [22]. Let us briefly summarize the main features of the 
computed flow structure. Figure El shows the three-dimensional surface of constant density 
p = 5 X V^'^^gjcw?. The region of the interaction between the stream and circumdisk halo 
is shown enlarged in Fig. El [22, Figs. 7, 8]. The left panel of Fig. (21 displays contours of con- 
stant density and velocity vectors, while the right panel of Fig. (21 is the so-called texture: a 
visualization of the velocity field by means of numerous tracks of test particles. 

The results presented demonstrate that the interaction between the circumdisk halo and 
the stream possesses all the characteristic features of an oblique collision of two flows. The 
resulting structure of two shock waves with a tangential discontinuity between them is clearly 
visible in Fig. (21 The region of the shock interaction between the stream and halo has a complex 
shape. The parts of the halo far from the disk have low density, and the shock due to their 
interaction with the stream lies along the edge of the stream. The shock bends as the gas 
density in the halo increases, finally assuming a position along the edge of the disk. Note that 
the solution for the cool case has the same qualitative features as the solution for hot outer 
parts of the disk [27]: the interaction between the stream and the disk is collisionless, the region 
of enhanced energy release is due to the interaction between the gas in the circumdisk halo and 
in the stream and is located outside the disk, and the resulting system of shocks is extended 
and can be considered as a "hot line". 

It follows from the above general features of the flow pattern that, at the interaction zone, 
the halo gas and stream gas pass through the shocks corresponding to their own flows, are 
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Figure 3: Contours of constant density and velocity vectors (left panel) and a visualization of the 
velocity field (right panel) in the region of interaction between the stream and the circumdisk halo, in 
the system's equatorial plane. This figure was first published in [22]. 
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Figure 4: Contours of constant density in the section of the flow structure perpendicular to the 
equatorial plane and passing through the accretor and the region of the halo's maximum thickening 
((/> = 0.7). 

mixed, and then move along the tangential discontinuity between the two shocks. Subsequently, 
the disk itself, halo, and inter-component envelope are formed of precisely this matter. The 
jump in the gas parameters after the passage of the shock leads to a the increase of density 
and temperature in the region between the shocks and, consequently, to the appearance of 
a pressure gradient along the z axis, perpendicular to the system's plane of rotation. As a 
result, the gas begins to expand vertically, increasing the z component of the velocity, until the 
pressure gradient is balanced by the gravitational force. 

The vertical gas pressure due to the presence of the z component of the velocity, together 
with the motion of the gas along the tangential discontinuity at the disk's outer edge, lead to 
a gradual increase of the thickening of the circumdisk halo (along the z axis). This thickening 
of the halo along the outer edge of the accretion disk is clearly visible in the three-dimensional 
constant-density surface shown in Fig. |21 as well as in Fig. |3J which displays const ant density 
curves in the section of the flow structure that is perpendicular to the equatorial plane and 
drawn through the accretor and the maximum thickening of the halo {(f) = 0.7). An additional 
illustration of the halo thickening above the disk is provided by Fig. El which shows fragments 
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Figure 5: Fragments of streamlines emerging from the neighborhood of Li. 



of computed streamlines emerging from the vicinity of Li. We can see that the streamhnes 
diverge after acquiring a vertical acceleration in the region between the shocks, then converge 
again toward the system's equatorial plane. The region of vertical acceleration is restricted to 
the hot-line zone. After passing this region, the gas has a large vertical velocity component 
that makes it climb higher, until its store of kinetic energy is exhausted. The point where the 
upward motion ceases corresponds to the maximum height, which is reached at phase ~ 0.7, 
i.e., already substantially outside the hot line. 

To quantitatively analyze the thickening of the halo above the disk, let us consider the 
behavior of the streamlines and the distributions of the gas parameters along them. In a 
cylindrical coordinate system with its origin coincident with the accretor {x = A,y = 0.0, z = 
0.0) and the angle measured from the point Li opposite to the direction of rotation of the 
matter (coincident with the system's direction of rotation), each point of the streamline is 
described by the coordinates {r,(j),z). The z{(f)) relations for four streamlines originating at 
points in the neighborhood of Li are presented in Fig. El We can see that, after entering 
the hot-line region (phase ~ 0.975), the streamlines begin to climb due to the increase of the 
velocity's z component. The phase relation of the vertical velocity for the same four streamlines 
is displayed in Fig. [7| When the gas emerges from the hot line (phase ~ 0.8), the force due 
to the pressure jump behind the shocks disappears, and the gas simply moves in the accretor's 
gravitational field. The vertical velocity of the gas becomes zero at phase ~ 0.7, corresponding 
to the maximum ascent of the streamlines. The height of the halo at this position reaches 
~ OMA. 

In addition to the primary maximum, the system also exhibits a minimum height at phase 




Figure 6: z((/>) relation for four streamlines. The streamlines originate in the vicinity of Li and have 
the coordinates (x^j , 0, zq). The zq values for each streamline are indicated in the plot. The shaded 
area corresponds to the region where Vz increases. 




Figure 7: Vz (</>) relation along the same four streamlines as in Fig. |S1 
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Figure 8: The density relation p((j)) along the same four streamlines as in Fig. El 

~ 0.3, when z ~ (0.005 -r- 0.006)A, and a secondary maximum at phase ~ 0.1, when the height 
of the halo is ~ O.OIA. Due to the effects of viscous dissipation, the minimum and secondary 
maximum are considerably less pronounced than the primary maximum. Eventually, viscosity 
leads to complete damping of the vertical velocity oscillations, so that the gas no longer possesses 
a significant vertical velocity when it next approaches the region of the interaction with the 
stream. 

Figure |H1 is the phase dependence of the density for the same four streamlines as in Fig. El 
The densities correspond to the mass-exchange rate in the RW Tri system {10~^MQ/yr [28]). 
We can see from these results that a strong density decrease is observed in the region of 
maximum ascent of the streamlines. This behavior of the density can lead to a displacement 
of the light-curve dips towards regions of higher p; i.e., towards higher phases. This should 
be manifest most strongly in short-wavelength observations, since higher densities are required 
to absorb harder radiation. It is also interesting that the gas density is much higher near the 
secondary maximum, at phases ~ 0.1, than near the first thickening. 

To evaluate whether the central source can be eclipsed by the thickening of the halo, we 
computed the optical depth of the halo matter, r. The optical depth is the product of the 
density, the layer's geometrical depth, and the absorption coefficient: r = pin. To estimate r 
near the region of thickening, let us use the densities presented in Fig. |H1 taking the charac- 
teristic linear size of the thickening to be its halfheight (Fig. We estimate the absorption 
coefficient using the approximation formulas from [29] with a temperature of ~ 13600°i^r and 
a characteristic density corresponding to the value at the thickening's half maximum. Our 
estimates of r are presented in Fig. O which shows the relations between the optical depth and 
the phase, t(0), along the same four streamlines as in Fig. (HI An analysis of these estimates of 
T shows that the thickening is optically thick, even at its greatest heights, and the value of r at 
the half maximum can reach ~ 10^. This means that the halo thickening is capable of eclipsing 
the central source, giving rise to the dips observed in the light curves of semi-detached binaries 
with stationary disks. 

Obviously, in stationary systems, when there is no collision between the stream and the 
outer edge of the disk, there can be no ricochet and no transfer of a substantial fraction of the 
stream's matter into the inner parts of the disk. However, the momentum of the stream matter 
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Figure 9: The optical depth relation, t(0), near the thickening (0 = 0.55 -^0.85) along the same four 
streamlines as in Fig. |EJ 
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is lower than that of the halo gas, so that some of the stream gas at great heights can slide 
down into inner orbits, appearing as if the stream is flowing around the outer edge of the disk. 
This sliding-down can take place only in the hotline region, where the momenta of the matter 
in the stream and in the halo have not yet become equalized. Indeed, when considering the 
streamlines that reach high heights from the stream's center fFigllUj). we can see that some of 
them leave the interaction region of the stream and halo, envelop the disk from the upper side, 
and get into its inner regions. The results of our computations demonstrate that the total flux 
of matter associated with this effect is small, and apparently cannot significantly influence the 
visibility of the central object. 

4 Conclusions 

Observations of semi-detached binaries with stationary accretion disks reveal dips in their 
light curves at phase ~ 0.7. These dips are usually explained by invoking the presence of a 
thickening at the outer edge of the accretion disk, without explaining how matter appears at 
considerable heights above the accretion disk in the case of a stationary interaction between 
the stream and disk. It is important to resolve the problem of the formation of this thickening 
of the halo above a stationary disk if we wish to correctly interpret the observations and more 
fully understand the nature of matter-exchange processes in these stars. 

The results of our three-dimensional numerical modeling of the matter-flow structure in 
semidetached binaries with stationary disks confirm our earlier conclusions [21, 22, 27] that the 
interaction between the stream and disk is coUisionless, the region of increased energy release 
is due to the interaction between the gas in the circumdisk halo and the stream and is located 
outside the disk, and the system of shocks that forms is extended and can be considered as 
a "hot line." The interaction between the circumdisk halo and the stream possesses all the 
characteristic features of an oblique collision of two flows that results in the formation of a 
structure consisting of two shock waves and a tangential discontinuity. 

In the hot-line region, the halo and stream gases pass through the shocks corresponding to 
their own flows, are mixed, and then move along the tangential discontinuity between the two 
shocks. During the interaction between the stream and the circumdisk halo, a considerable 
fraction of the matter acquires a vertical acceleration. The vertical motion of the gas due to 
the z component of the velocity, together with its motion along the tangential discontinuity 
at the outer edge of the disk, results in a gradual growth of the thickness of the circumdisk 
halo. The region of vertical acceleration is restricted to the hot-line zone, and its angular size 
does not exceed ~ 65°. However, once it has passed this region, the gas has a sufficiently 
large vertical velocity component to rise until its store of kinetic energy is exhausted. The 
point where the upward motion ceases corresponds to the maximum height, which is reached 
at phase ~ 0.7, already substantially outside the hot line. The thickening extends appreciably 
higher than the scale height of the disk, reaching values ~ 0.04A (this corresponds to a ratio 
of the thickening's height to the distance to the accretor > 0.1), and its angular size exceeds 
~ 130°. Our computations also show that, in addition to the primary maximum, the system 
also displays a height minimum at phase ~ 0.3, when z ~ (0.005 0.006)^4, and a secondary 
maximum at phase ~ 0.1, when the height of the halo is ~ O.OIA. 

Our analysis of these results leads us to conclude that the dips in the light curves of semi- 
detached binaries with stationary disks (i.e., in the absence of a collisional interaction between 
the stream and disk), can be explained by the formation of a thickening of the halo above the 
outer edge of the disk. The origin of this thickening is described well by the hot-line model, and 
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its quantitative features arc consistent with observations. The proposed model can be applied 
to both semi-detached systems (LMXBs, cataclysmic variables) in their stationary state and to 
dwarf novae in outburst, provided that the outer parts of the disk are not strongly distorted in 
the outburst. 
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